DMFT vs Second Order Perturbation Theory for the Trapped 2D 
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In recent literature on trapped ultracold atomic gases, calculations for 2D-systems are often done 
within the Dynamical Mean Field Theory (DMFT) approximation. In this paper, we compare 
DMFT to a fully two-dimensional, self-consistent second order perturbation theory for weak in- 
teractions in a repulsive Fermi-Hubbard model. We investigate the role of quantum and of spatial 
fluctuations when the system is in the antiferromagnetic phase, and find that, while quantum fluctu- 
ations decrease the order parameter and critical temperatures drastically, spatial fluctuations only 
play a noticeable role when the system undergoes a phase transition, or at phase boundaries in 
the trap. We conclude from this that DMFT is a good approximation for the antiferromagnetic 
Fermi-Hubbard model for experimentally relevant system sizes. 

PACS numbers: 67.85.-d, 75.50.Ee, 71.10.Fd 
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INTRODUCTION 

The study of ultracold, trapped atomic gases on a lat- 
tice as an emulation of tight-binding Hamiltonians has 
been an active field of research for the past 15 years. Af- 
ter the realization of BECs of alkali atoms in 1995 [H, 
many theoretically predicted quantum phenomena, such 
as the bosonic superfiuid-Mott insulator transition[3| and 
supercxchangcQ, have been observed experimentally. 
The recent achievement of a fermionic Mott insulator in 
a Hubbard Hamiltonian]^ Q has steered the focus to- 
wards understanding the fermionic Hubbard model and 
its phase diagram, which may have a relevance for high 
temperature superconductivity^,!!]. An important chal- 
lenge is to reach the Neel temperature, in order to realize 
antiferromagnetic phases. 

It is common practice, when doing calculations on 
two-dimensional trapped atomic gases, to use DMFt[|, 
which includes only spatially local fluctuations. While 
DMFT is known to be a good approximation in three 
dimensional problems [lol [ll|. it is expected to perform 
much more poorly in two dimensions, since non-local fiuc- 
tuations gain importance [l^, [13 • 

We have therefore implemented a self-consistent sec- 
ond order perturbation theory for the Fermi-Hubbard 
Hamiltonian. The self energy expansion, which is de- 
scribed further on, includes local as well as non-local 
fiuctuations of the system. This gives us the tools to 
compare the effects of quantum and spatial fiuctuations 
for the model, up to second order in the interaction. 

A two pscudospin-species Fermi mixture in a magneto- 
optical trap realizes a Hubbard model with broken lattice 
symmetry, 



where (ij) denotes the sum over next neighbors, and 
C;^ are the creation and annihilation operators for the 
respective lattice point and spin and h^^ — c\^c.^^ is the 
number operator. Also, t is the next-neighbor hopping 
amplitude, and U the on-site interaction, which is cho- 
sen to be repulsive, U > 0, in this paper. The term 
Vi^ describes the harmonical trapping potential which 
breaks the lattice translational symmetry, and the spin- 
dependent chemical potential is ^^r- The energy scale for 
this work is the hopping amplitude, t {t = 1). 

It is necessary to use approximative schemes to solve 
this Hamiltonian. In this paper, three different ap- 
proximations are compared, namely Hartree-Fock the- 
ory, DMFT, and self-consistent perturbation theory up 
to second order in the interactions. It is known from 
Hartree-Fock theory fl3| that, for the case of an imbal- 
ance, 7^ /^j,, the ground state is a canted antiferromag- 
net, with a C/(l)-symmetry for all rotations around the 
z-axis. For the case of balanced spin species, = /^|, the 
system is symmetric under all rotations of SU(2). While 
the code used can easily reproduce both balanced and 
imbalanced cases, this paper will be constrained to bal- 
nced systems, for simplicity. These are sufficient for an 
identification of key effects of second order diagrams. 

This paper is organized as follows. First, we will give 
a short introduction to the approximation schemes used 
here. Then, we will discuss results for the trapped sys- 
tem, focusing on the differences of the approaches. The 
role of non-local diagrams is then investigated further, 
and finally a short survey of finite size effects, which are 
noticeable in our calculations, is presented. 



THE DIFFERENT APPROXIMATIVE SCHEMES 

In Hartree-Fock theory (HF), quantum fluctuations 
are completely neglected. This is expressed through the 



2 



transformation 
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This is equivalent to a self-consistent perturbation 
scheme of first order [c.f. Figs. [T] a) and b)], with a 
local and frequency-independent self energy (SE), which 
is therefore diagonal in real space representations: 



(3) 



The neglect of quantum fiuctuations leads to a systemati- 
cal exaggeration of both critical temperature and ground 
state staggered magnetization in all dimensions. 
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FIG. 1: First and second order SE diagrams, a) Hartree-term 
b) Fock-term c) density-density term d) exchange term 

Dynamical Mean Field Theory (DMFT) tries to 
correct for this shortcoming by including some quantum 
fluctuations, namely all local ones. The resulting self 
energy continues to be diagonal in real space, but now 
has a frequency dependence: 



HF 



(4) 



DMFT is usually viewed as the limit of infinite dimen- 
sions, and as such is a useful method for high [d > 3) 
dimensions. It is generally expected, however, that in 
lower dimensions the usefulness of DMFT is strongly re- 
duced, because nonlocal processes gain in weight. DMFT 
is a priori a nonperturbative expansion: the SE includes 
diagrams of all orders of perturbation. For the sake of 
our comparison, we have only included contributions to 
the local SE up to second order. We expect this to deviate 
only slightly from the full DMFT, for the same reasons 
given below for full second order calculations. In the dia- 
grammatic language, this corresponds to an inclusion of 
diagrams a) to d) in Fig. [1] but only those whose Green 
functions' lattice indices are equal, Gia-,iK for all Green 
functions in the diagram. 

Finally, the self-consistent second order pertur- 
bation theory (20PT) includes local and nonlocal dia- 
grams. The SE is now nondiagonal in real space. 



^20PT 

i(7,jK, 



(5) 



The numerical effort grows exponentially with every or- 
der of interactions in the perturbativc series, making 



it numerically difiicult to include more than second or- 
der diagrams. Because we stay in the weak interaction 
regime, and higher orders of interactions are suppressed 
by at least U^/W'^ - where W = At is the typical half- 
bandwidth of the two-dimensional Hubbard model - sec- 
ond order calculations should include all dominant quan- 
tum fluctuations and, therefore, return quantitatively 
correct results. 



THE INHOMOGENEOUS SYSTEM 



a) density distribution, Hartree-Foci<, U=3.5, V=0.056 



b) magnetization, Hartree-Focl<, U=2.5, V=0.056 
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c) density distribution, DIVIFT, U-3.5, V-0.056 
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d) magnetization, DMFT, U-3.5, V=0.056 
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density distribution. 20PT, U=3.5, V-0.056 
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f) magnetization, 20PT, U=3.5, V-0.056 
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I) density distribution, p-50, U-2.5, V-0.0034 h) magnetization, p=50, U=2.5, V=0.0034 
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FIG. 2: (Color online) Cuts through the lattice for differ- 
ent approximations, displaying the occupation and the stag- 
gered z-magnetization. Figs, a) to f) show almost the entire 
occupied region, with trap potential V = 0.56, while Figs, 
g) and h) show the center of a system, with V — 0.0017. 
The HF results show significantly higher critical temperatures 
and higher ground state staggered magnetizations, as well as 
sharper phase transitions, than the other two approximations, 
20PT and DMFT, which differ only slightly, although con- 
sistently, from another. 



For the Fcrmi-Hubbard model on a square lattice, the 
additional trapping potential leads to a coexistence of 
different phases in the system, depending on the local 
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density: For local densities lower and significantly higher 
than half filling, a paramagnetic phase is present. In 
regions close to half filling, a Neel state is present. While 
at high temperatures the Neel state adapts to the filling 
as given by the trap potential, when temperatures are 
lowered the Neel phase enforces half filling over a broad 
range, reminiscent of the occurrence of phase separation 
in a doped homogeneous system[l5|. Unlike the Mott- 
plateau, which is a fluctuation effect not reproducible in 
HF-calculations [l^l, the Neel-plateau can be seen even 
at HF-level. 

A domain wall boundary appears as soon as the sys- 
tem crosses into the Neel-plateau region IJ]. At the 
new boundary, the antiferromagnetic order parameter 
{—ymz{i) changes its sign. The boundary can also be 
seen in the occupation number as a ring of higher occu- 
pation extending from the region of half filling [Figs. [H 
g), h)]. Notice that the higher occupation is suppressed 
again towards half filling by the second antiferromagnetic 
domain, working in opposition to the trapping potential. 
This new Neel domain requires the occupied part of the 
lattice to extend to a radius of at least 15 to 20 lattice 
points around the trap center, otherwise the lattice res- 
olution is too small for the new domain to exist. 

Our calculations show that quantum fluctuations in- 
crease the inverse Neel temperature, /3jv, while decreasing 
the Neel order parameter, c.f. Fig. [2] In HF, P^^{U = 
3.5) is 1.5, while P§^'^^{U = 3.5) and Pjy^^{U = 3.5) 
are 3. The (effectively) ground-state staggered magneti- 
zation at [/ = 7 is reduced from A"^{/3 = 30) = 0.319 to 
^DMFT (30) = 0.235 and A2O^'^(30) = 0.226, by 30%. 
The progression from low to high temperatures can be 
seen in Figs. [2]a) to f). In HF, the magnetization remains 
sizeable up to comparatively high temperatures, and the 
crossover to the low-T Neel state occurs at (5^^ = 4 [c.f. 
Fig. [2] a) , where the plateau disappears for high-T cal- 
culations]. In Figs. Elc) to f), the results for DMFT and 
20PT arc shown. The crossover behavior is smeared out 
over a wide range of temperatures, starting to set in at 
/3 = 8 and being fully developed only around /3 = 18. 



THE ROLE OF NON-LOCAL DIAGRAMS 

In this section, we want to concentrate further on the 
conformance of 20PT and DMFT results. To this end, 
we take a look at finite systems with periodic boundary 
conditions, leaving out the trapping potential, = 0, for 
translational invariance. This gives access to the effects 
of local and non-local diagrams without parallel effects 
from trap curvature or boundary potentials. 

Fig. [3] offers a typical progression of 20PT corrections 
to the systems: At very weak interactions, no order is vis- 
ible. The critical interaction U^^{13) is reached in HF 
before U^^^"^ {P) is reached in DMFT, which in turn is 
closely followed by U]p^'^ {j3). In medium- low interac- 
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FIG. 3: (olor online) Comparison of HF to 20PT and DMFT 
calculations on the [/-axis. The system is a 16 by 16 lattice, 
with periodic boundaries at half filling and low temperatures 
(/3 = 50). 20PT and DMFT difler only marginally. 



tion ranges, the suppression reaches a maximum, usually 
around 2 < [/ < 3, and then slowly weakens at higher 
interactions. This last piece of information must, how- 
ever, be treated with caution, as our method is of limited 
validity at higher interactions. 




FIG. 4: (Color online) Ratio of local to non-local elements in 
the self energy field. Fig. a) shows the ratio for two different 
temperatures T\ = 0.02, which is strongly antiferromagnetic, 
and T2 = 0.11, which is only very weakly magnetized. The 
weight of non-local elements never exceeds about 20%, and 
sinks to about 5% in an area around the Fermi energy. Fig. 
b) shows a selection of the weight of diflerent shells in real 
space around purely local elements; Shell 1 thus means local 
and next-neighbor self energy elements, shell 5 all elements 
with exchange between lattice site up to 5 jumps (Manhattan 
metric) from the original site. 

Again in Fig. [3J the close proximity of 20PT and 
DMFT results when sampling systems with different in- 
teraction strength U at constant temperature and filling 
can be seen. While DMFT is shghtly above 20PT for 
the entire graph, the difference is only 1-4% at interac- 
tions from U = 1.25 onwards. The only clear deviation is 
near the critical interaction tV, where 20PT shifts the 
critical interaction further from U^^{f3) than DMFT. 

The close agreement of the data can also be seen in the 
self energy. To ascertain which role is played by non-local 
fluctuations, we compared the weight of local and non- 
local self energy diagrams in a 20PT calculation. The 
calculation was done by adding the moduli of all local self 
energy matrix elements, and normalizing to the summed 
moduli of all self energy diagrams for each frequency. The 
result, shown in Fig. 21 is that off-diagonal elements con- 
tribute only weakly to the antifcrromagnct. Especially 
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around the Fermi frequency, the local diagrams make up 
95% of the self energy's weight. The picture is essentially 
the same for high and low temperatures, and for different 
relevant interaction strengths. 

We conclude from these graphs that, for antiferromag- 
netic systems, the influence of non-local fluctuations is 
weaker than generally expected, and that DMFT is a 
valid approximation in the trapped weak coupling 2D- 
antiferromagnet. The only notable divergences are near 
phase boundaries, both spatial boundaries, as in Figs. [2] 
g) and h), or in parameter space, as in Fig. |31 

We were able to investigate lattice sizes of up to 18 
by 18 atoms. Since this is smaller than common experi- 
mental extents, we include an investigation on finite size 
effects on the system in this paper. For these runs, we 
again chose periodic boundary conditions and no har- 
monic potential, in order to investigate only effects com- 
ing from the increasing number of nonlocal diagrams and 
lattice points, and eliminate effects stemming from differ- 
ent curvature of the trapping potential in the discretized 
lattice, and from the constant boundary condition of the 
experimental setup. 





FIG. 5: (Color online) a) The relative suppression of the an- 
tiferromagnetic phase through second order diagrams as com- 
pared to HF results. The labels refer to the system size, and 
all calculations are done at /? = 50 and half filling in a periodic 
system. The larger the system, the stronger the suppression, 
b) The deviation of DMFT and 20PT results for exemplary 
system sizes. Shown are the regions where both approxima- 
tions have appreciable staggered magnetization. For all sys- 
tem sizes investigated, no noticeable trend can be made out. 

The effect of finite system sizes is pronounced, as can 
be seen in Fig. [5] a): The critical interaction Um{P) is 
shifted by 20PT; the larger the system the stronger the 
shift. The amplitude of the suppression also increases 
with the system size, so that larger systems reduce an- 
tiferromagnetic order, as expected from the Mermin- 
Wagner theorem. As a direct consequence of this, we 
expect experimental results to show a stronger reduction 
of antiferromagnetic order, and thus a lower critical tem- 
perature, than what we have calculated here. 

In Fig. [S] b) , we show the relative deviation of DMFT 
to 20PT runs for different system sizes. No trend in lat- 
tice extents can be made out. We deduce, that the spatial 
fluctuations converge quickly for increasing shells around 
the local approximation, and do not change appreciably 
for system sizes of more than 10 by 10. All flnite size 
effects in Fig. [5] a) thus stem from local fluctuations. 
CONCLUSIONS 



ond order self-consistent perturbation theory, and com- 
pared these results to HF and DMFT approaches. 



We find that the inclusion of quantum fluctuations to 
HF calculations strongly shifts critical temperatures and 
order parameters to lower values, as expected. Inclusion 
of non-local parts to the self energy changes little over 
a large range of values. Only near critical points and 
boundaries do noticeable variations appear between sec- 
ond order perturbation theory and dynamical mean field 
theory. This is refiected in the small weight of non-local 
elements in the self energy. 



We then looked at the progression of the suppression 
with increasing lattice size. The suppression through 
quantum fiuctuations increases monotonically with the 
size, showing pronounced finite size effects. Spatial fluc- 
tuations, on the other hand, quickly converge at small 
system sizes, and are fully represented in the systems we 
investigate. 
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In summary, we have investigated the antiferromag- 
netic phase of ultracold atoms trapped on a lattice in sec- 



